Variable-frequency-controlled coupling in charge qubit circuits: Effects of microwave 

field on qubit-state readout 
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To implement quantum information processing, microwave fields are often used to manipulate 
superconuducting qubits. We study how the coupling between superconducting charge qubits can 
be controlled by variable-frequency magnetic fields. We also study the effects of the microwave fields 
on the readout of the charge-qubit states. The measurement of the charge-qubit states can be used 
to demonstrate the statistical properties of photons. 



I. INTRODUCTION 

Superconducting quantum circuits are good candidates 
for implementing quantum information processing [l|, Q • 
To construct universal quantum computing, controllable 
couplings between any pair of qubits are required. Theo- 
retical methods for switchable couplings in charge-qubit 
circuits have been proposed by changingthe amplitude 
of the bias magnetic flux, e.g., in Refs. [JM II] ■ However, 
in experiments, it is much easier to produce fast and pre- 
cise frequency shifts of the radio-frequency (RF) control 
signals, as opposed to changing the amplitude of the dc 
signal. Methods using variable-frequency-controlled cou- 
plings in superconducting flux-qubit circuits have been 
studied [5[ theoretically, comparing with the coupling ap- 
proach using the dressed states [& 0, Hi ■ In this scheme, 
the two qubits can be coupled to (or decoupled from) 
each other by modulating the frequencies |f| of exter- 
nally applied variable-frequency magnetic fields to match 
(or mismatch) the combination of frequencies of the two 
qubits. The coherent oscillations and conditional gate 
operations of two superconducting charge qubits with 
always-on coupling have been demonstrated [a] exper- 
imentally. Therefore, the next step for charge qubits 
would be to design superconducting quantum circuits 
with switchable couplings. 

Here, we first generalize our approach [j| using the 
variable-frequency-controlled coupling in flux qubit cir- 
cuits to the charge qubit circuit proposed in Ref. 
This proposal has the following advantages: (i) the cou- 
pling between different charge qubits can be implemented 
fast by changing the frequency of the externally applied 
classical field; (ii) these proposed charge qubits always 
work at their optimal points, and thus the qubits are 
mostly immune from charge noise [lol ]. produced by un- 
controllable charge fluctuations; (iii) no additional circuit 
is needed to realize this controllable coupling. 

Besides the controllable coupling, measuring the qubit 
state is also a very important step in quantum informa- 
tion processing. In superconducting quantum circuits, 
microwave fields are often used to implement quantum 



information processing. Here we focus on how microwave 
fields affect the readouts of the qubit states. In particu- 
lar, we explore the effect of quantized fields with differ- 
ent statistical properties on measurement results of the 
qubit states when the charge qubits are placed inside a 
microcavity, e.g. a three-dimensional cavity (ill . Il2| or a 
superconducting transmission line (l3j . 

The paper is organized as follows: in Sec. II, we 
generalize the variable-frequency-controlled coupling ap- 
proach in flux-qubit circuits [5] to that in charge-qubit 
circuits [4]. In Sec. Ill, we study the effect of the clas- 
sical and quantized microwave fields on the readout of 
the qubit states. In Sec. IV, we compare the classical 
and quantum treatment of the large Josephson junction. 
Finally, conclusions are presented in Sec. V. 



II. HAMILTONIAN WITH 
VARIABLE-FREQUENCY CONTROLLED 
COUPLINGS 

We first verybriefly review the model Hamiltonian pro- 
posed in Ref. [4j for two coupled superconducting charge 
qubits by sharing a large Josephson junction (JJ) (see 
Fig. 1). The large J J is classically treated and its charge 
energy E c0 is neglected. The dc biased magnetic field $ e 
is externally applied through the area between the large 
J J and the first qubit. Each qubit is also biased by a dc 
voltage Vxi via the gate capacitance Ci (i = 1, 2). The 
Hamiltonian of the superconducting circuit is [J] 
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with Ei(Vxi) = E c% {n t - C l V Xl /2e) 2 . Here E cl = 
2e 2 /(Ci + 2Cji) and Eji are the charge and Josephson 
energies of the ith charge qubit. Ejq is the Josephson 
energy of the large J J. The number rii of excess Cooper 
pairs in the superconducting island is canonically conju- 
gate to the average phase drop (fii = (ifiA + of the 
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FIG. 1: (Color online) Schematic diagram of two charge 
qubits coupled by a (left) large Josephson junction (JJ) with 
coupling energy Ejo and capacitance Cjo- For the ith charge 
qubit (where i = 1, 2), a superconducting island (denoted by a 
filled circle) is connected to two identical small JJs (each with 
coupling energy Eji and capacitance Cji). Also, this island 
is biased by the voltage V% = Vxi + Vg (i) via a gate capaci- 
tance d, where Vxi is a static (dc) gate voltage and V g (t) 
is a time-dependent (ac) microwave gate voltage. Moreover, 
a static (dc) magnetic flux <I? e plus a microwave-field-induced 
magnetic flux &(t) (ac) are applied to the (yellow) region be- 
tween the large JJ and the first charge qubit. 



to the superconduction island of the ith qubit via its gate 
capacitance, and an additional variable-frequency (ac) 
magnetic flux <!> (t) = <I> C sin (cut) is also applied through 
the area between the large JJ and the first charge qubit 
(see Fig. [T]). To make our proposed charge-qubit more 
immune from the uncontrollable charge fluctuations, it is 
also assumed that two charge qubits work at their opti- 
mal points, i.e. the applied dc voltages Vxi satisfy the 
condition £i(Vxi) = 0. Considering these conditions, the 
Hamiltonian in Eq. ((2]) becomes 
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where e { Q t} = E cl (C t V gt /2e) and g l = 2E, n sin (7r$ e /$ ) f ■ 
The parameters gvi and £ are given by 
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ith charge qubit. The phase drop across the large J J is 
7. Considering that the critical current Iq = 2itEjq/Qq 
of the large JJ is much larger than the critical currents 
Id = 2nEji/$ of the charge qubits, the phase 7 across 
the large JJ is very small. We can expand the func- 
tions of the phase drop 7 in Eq. (fT]) into a series and and 
retain the terms up to second order in the parameters 
Vi — {Id/Io) < 1. In this case, Eq. (fl]) can be reduced to 
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Here cp c = 2tt$ c /3> an d Jn is the nth-order Bessel func- 
tion of the first kind. 

In the rotating reference frame at the frequency LJ g i 
about a x ^ , the Hamiltonian in Eq. ((4]) is rewritten as 



in the spin-i representation based on the charge states 
|0) = I t) an d |1) = ||) that correspond to zero and 
one excess Cooper pairs in each Cooper-pair box, where 
ei(Vxi) = ±E ci (CiV Xi /e-l), and 
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with i, j = 1, 2 (i ^ j). The coupling constant X12 
between the two charge qubits is 



Xi2=LjI c iI C 2 sin 
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To eliminate the a x term in Eq. ([5]), the frequency 
of the microwave field applied to the gate capacitance 
is set as hu> g i ~ Eji. Furthermore, we can tune the 
flux $ e so that the coupling strength \\2 is less than the 
coupling strength gi2- Also, we tune the gate voltage V g i 

so that the large detuning condition — Eq 1 ' = A ^> 

X12 can be satisfied. Under this condition, the always-on 
interaction X12 is negligibly small, and the Hamiltonian 
in Eq. © is reduced [lj to 
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where Lj = <J> /27rio is the Josephson inductance of the 
large JJ. It is clear that the coupling between the two 
qubits is realized via this effective inductance. 

Now, we study how to apply our variable-frequency- 
controlled approach (5j to the above charge-qubit cir- 
cuits 0. We assume that besides the dc voltages Vx% 
and the dc magnetic flux $ e , an ac microwave voltage 

Vg l \t) = Vgi cos (ujgit) with the frequency u> g i is applied with hio\ = £q — \ , hbj 2 — s [ Q z> + x' , and x' = X12/2A 
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Let us discuss how the interaction between two qubits 
can be switched on and off via Eq. by changing the fre- 
quency u> of the variable-frequency magnetic flux <fr(t) = 
& c s'm(u}t). Equation © shows that the two qubits are 
approximately decoupled from each other when there 
is no applied ac magnetic flux $(£). However if the 
frequency uj of <I>(i) is tuned to satisfy the condition 
u = ui\ + 0J2, then two qubits can be simultaneously 
flipped by the variable-frequency magnetic flux via the 
interaction Hamiltonian 



t r (1) (2) 
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where the contributions of other fast oscillating terms are 
negligibly small. If the frequency uj of $(t) satisfies the 
condition uj — u)2 — u>\ , then one qubit can be flipped by 
another with the help of the variable-frequency magnetic 
flux through the interaction Hamiltonian 
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after neglecting other fast oscillating terms. 

A single-qubit operation can also be implemented via 
the variable-frequency magnetic flux For example, 

if uj = uii, or ui = u>2, then the first or second qubit can 
be selectively rotated around the x axis. When there is 
no variable-frequency magnetic flux, a rotation around 
the z-axis can be implemented for each qubit. Therefore, 
any logic gate (see, e.g., Ref. [THh can be realized by using 
single-qubit operations and two-qubit operation via the 
Hamiltonians in Eqs. (J7J and ([5]). 



III. EFFECT OF MICROWAVE FIELDS ON 
SUPERCURRENTS IN THE MEASUREMENT 
OF QUBIT STATES 

Above, we have shown that the interaction between the 
two qubits can be switched on and off using a variable- 
frequency magnetic flux. Two-qubit operations can be 
implemented, and entangled states between two qubits 
can also be generated, using Eq. ([7J or {5J. To imple- 
ment the readout of two-qubit states, we need to calcu- 
late the circulating supercurrent I contributed by the two 
qubits [ij . The operator of the supercurrent I of the two 
qubits is given by 



/ = sin (^] (/,,,,!"- Ua) 1 
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For any given state (e.g., of two qubits, the super- 
current can be obtained by 

7=<*|/>>. (10) 

Note that two-qubit operations are always related to the 
microwave fields. The supercurrent / might be differ- 
ent for different microwave fields with different statistical 
properties. Below, we study how the different microwave 
fields affect the supercurrent /. 



A. Classical microwave field 

We now focus on two-qubit entangled states, created 
from the ground state 1 51,172) via the two-qubit interac- 
tion Hamiltonian in Eq. ([7]) . For these created entangled 
two-qubit states, the contribution of the average values 
of single-qubit operators (i = 1, 2) to the supercur- 
rent is zero, and the supercurrent I is only determined 
by the two-qubit operator ai^a^ as follows 



{I) 



■ sm 



27T$ e 
$0 



1 



(11) 



Here, for simplicity, the two qubits are supposed to have 
identical critical supercurrent Z c i = I c2 = J c , and then 
?]i = T)2 = ij- The quantum fluctuation of the total su- 
percurrent I is 



(Air = (n (if, 
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which can be further given by 
7] 2 2 . o ( 2tt$ 
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Considering that the ratio rj is small, the first term in 
Eq. (fl3|) can be neglected in the following calculations. 
In this case, the supercurrent fluctuation has a simi- 
lar behavior to the supercurrent (I) of Eq. (Ilip and 

(A J) 2 oc (I) , when entangled two-qubit states are created 
from the ground state 131,32) through the Hamiltonian 
in Eq. (J7J. For convenience, we define a reduced quan- 
tity k to describe the supercurrent (I) and supercurrent 
fluctuation (A/) 2 as 



(AJ)2 



2/2sin 2 (^T$ e /$o) 
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If two charge qubits are initially in an entangled state 
cos0|<7i, 32) +sin0e l ^|ei, 62) and the evolution of the two 
qubits is governed by the Hamiltonian in Eq. (J7J, then the 
reduced supercurrent or supercurrent fluctuation k(t) 
can be given by 

k c = 1 + sin(20) cos (r) cos 4> + 008(20) sin (r) , (15) 

which means that k c is an ac signal. Here, r = 
|<7i2|£, with the evolution time t. If the initial state is 
cos6\g\,g2)+sva.6\e\, e 2 ), then k c = l+sin(r+2#). When 



the evolution time To 

4. / (!) ( 2 ) 

gives rise to ( tri <ri 



-28 + 2nir 



0, which 



— 1. Thus, in this case both 

the total supercurrent and the supercurrent fluctuation 
become zero. 
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B. Quantized microwave field 

Now let us consider the case when the variable- 
frequency magnetic flux $ c cos(wi) is replaced by a quan- 
tized magnetic flux, <fr q a + + $*a, with frequency lu — 
lui+lu2- Following the same way as the above derivation 
of Eq. l[7]). we can obtain an interaction Hamiltonian Hj 
between the quantized magnetic flux and the two charge 
qubits 



C* C 1 ) ( 2 ) f1(f\ 

■£ u aayay (16) 



where 



2TT$ q LjI cl I c2 . {2w$ e 
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This model indicates that one photon can flip both qubits 
simultaneously. 

We now consider that the two qubits are initially in 
the state cosfl \g, g) + sin0 e l ^\e,e) and the quantum field 
is initially in a state D(n) \n), here D(n) will be given 
below for a given state. From the Hamiltonian (HHJ), the 
total system evolves to 
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+ cos 0£>(O) | <?,<?,()), (18) 

a„(r) = cos(r\/n+l ) sin 6 D(n) 

— sin(rvn + l ) cos 9 D(n + 1), 

frn( T ) = cos(tv / w + 1 ) cos 9 D(n + 1) 

+ sin(rVn + 1 ) sin 9 e # D(n), 



with the rescaled dimensionless time r = |^ 12 1 i- Using 
Eq. {Tip and Eq. (TT8"j) , at the time r, the reduced su- 
percurrent expectation value or supercurrent fluctuation 
Kq(r) in the case of the quantized field is 



k q (t) = 1+2 Re |wo(t) cos 0.0(G) - > [„„ ( r) r„ ( r )\ 
where 
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w„(t) = cos( 

+ sin(rVn+ 1) sinfl e i0 £>(n). 



Equation (119|) shows that the supercurrent ex- 
pectation value (7) consists of a dc component 
— (r)I c /2) sin(27r$ e /$o) and different ac components, 




FIG. 2: The photon number distribution P(n) with average 
photon number n — 7 of (a) coherent state \a); (b) superpo- 
sition of coherent states (\a) + \ — a))/N+; and (c) squeezed 
vacuum state |0, (). 



which are modulated by time-dependent factors, e.g. 
cos(rv / « + !)• 

We further specify that the quantized field is initially in 
several different quantum states [l6| , e.g., (i) the coherent 
state 



V2 E 



a 



(20) 



with a = ^/ne ltp ; (ii) the superposition of two coherent 
states 



N. 



-(\a) 



+ 



v 2n 



\2n) , (21) 



with N+ = V2(r 
state 



^/(2n)!cosh 
~ 2 "); and (iii) the squeezed vacuum 



|o,o = E 



n\\/coshr 



s l/3 tanh(r/2)]"|2n), (22) 



with the squeezing parameter £ = re 1 ^ and n = sinh 2 r. 
Here, n is the average photon number. The photon num- 
ber distributions P(n) (e.g., P(n) = \D{n)\ 2 = |(n|a)| 2 
for a coherent state) of the above three states are shown 
in Fig. [2 Physically, coherent states display Poisso- 
nian distribution and the fluctuations of both quadra- 
ture components are equal to the standard quantum fluc- 
tuation limit, AXi = AX 2 — 1/2. Squeezed states 
have sub-poisson distribution and the fluctuation for one 
of the quadrature components can be squeezed, e.g., 
AXi < 1/2. 
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If the qubits are initially in the ground state |<?i , 52) 
and the quantized field is initially in the coherent state 
I a), the reduced supercurrent expectation value or the 
supercurrent fluctuation is obtained from Eq. (II!: 



Ko(r) = 1 — 2cos</?e 



»=0 



A(n) sin(r\/n + 1) cos(Ti/n), 

(23) 

with A(n) = fi(™+3) /(nly/n + 1). Equations (JTSJ) and 
(|23|) show that the supercurrent expectation value (I) 



and the supercurrent nuctuation (A/)) 2 are very sensi- 
tive to the phase <p of the coherent state. If <p = ir/2, 
K q has only a dc component; however, when ip 7^ 7r/2, n q 
consists of many different ac components. 

If the qubits are initially in the ground state | (71,52), 
but the quantized fields are initially in the squeezed vac- 
uum states or superposition of coherent states, then from 
Eq. (fl"9|) , K q (r) is given by 

00 

K q (r) = 1 — 2 Re [B(n) sin(rV n + 1) cos(r-\/ra)] , 

n=Q 

(24) 

withB(n) = D*(n+l)D(n). Because of D(2n + 1) = 0, 
there are B{n) = and K q {r) = 1. Thus, the oscillatory 
evolution disappears. 

From Eq. (|14p . we know that the macroscopic super- 
current expectation value (I) can be described by n q . 
Figure ^ shows that the supercurrent of the charge 
qubits are different with the same initial qubit state 
(Iff, 5) + \ e T e ))/^ but with different initial states of the 
quantum field. From Eq. (fT9|) , in the case of the coherent 
state \a) with the phase ip = tt/2, the total supercur- 
rent (/) displays a sinusoidal-like evolution, as shown in 
Fig. [3]^ a). However, when tp = 0, the total supercurrent 
is shown in Fig. [3jb). If the quantized field is initially in 
a superposition of coherent states, the total supercurrent 
(I), as shown in Fig. [31(c), demonstrates the collapse and 
partial-revival phenomena. In the case of the squeezed 
vacuum state, the total supercurrent approximately dis- 
plays an ac current with a quasi-periodic evolution, which 
is demonstrated by Fig. [21(d). All irregular oscillations 
of the supercurrent expectation or supercurrent fluctu- 
ation reflect the coherent interference that comes from 
the coherent superpositions of the different photon num- 
ber states. The different initial photon states result in 
different output of the measurement of the charge-qubit 
states. Therefore, the measurement of the charge-qubit 
states can demonstrate the statistical properties of the 
photons, and charge qubits could be served as photon 
detectors. 



IV. QUANTIZATION TREATMENT ON LARGE 
JOSEPHSON JUNCTION 

In the Hamiltonian (fTJ) , the Josephson energy term 
E c qN 2 of the large JJ is neglected and the large JJ acts 



* 1 



y 1 



* 1 - 



y 1- 





1 I 1 I 1 


1 1 1 1 


1 1 1 1 1 
; (c) 

-k — ^§11^ 

II 


1 1 1 1 

aaaAAAA/WVWV v aaam 


' 1 ' 1 

1 1 1 1 


■ 1 1 1 1 

/VVW 

1 1 1 1 1 



10 



20 



30 



40 



50 



FIG. 3: Evolution K q (r) of the reduced total supercurrent 
expectation and reduced supercurrent fluctuation from the 
initial qubit state (\g,g) + \e,e))/y/2 in the presence of the 
quantum field initially in (a) coherent state \a) with the phase 
tp = 7r/2; (b) coherent state ip = 0; (c) superposition of coher- 
ent states (|q) + | — a))/N+ with p — 0; (d) squeezed vacuum 
state |0,£}- The irregularity of oscillations originates from 
the interference effect of the photon component of the above 
states. 



as an effective inductance Lj 0, [13] ■ We now consider 
a quantum mechanical treatment for the large J J. Con- 
sidering the additional term of charge energy E c oN 2 , the 
Hamiltonian of the large JJ can be written as 



H = E c0 N- - Ejq cos 7, 



(25) 



with the charge energy E c q and the excess Cooper pairs 
N. Because the large JJ works in the phase regime, the 
spectrum of the large JJ is approximately equivalent to 
a harmonic oscillator Hq = hiOpd^a, with the plasma fre- 
quency 



The bosonic operators a and arc defined by 

1 



and the phase drop 7 is expressed as 

1 



7 = -(« f + a) , 



(26) 



(27) 



(28) 



with c = (_Bj°' ) /2i?c ' ) ) 1 ^ 4 . Due to the large critical super- 
current of the large J J, one can expand the phase drop 7 
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in Eq. (QJ into a series and retain terms to the first order 
of 7. Finally, a spin-boson interaction between the two 
charge qubits and the large JJ is achieved: 

H = £ [e^Jo-W - E, H cos j of 
2 

+^0+0 + [giocrf (a + + a) , (29) 



where gio = — {Ej%/2<;) sin (7r$ e /<I>o). We assume that 
the plasma frequency w p of the large JJ is much larger 
than the splitting of the qubits. Thus the large JJ is al- 
ways in the ground state when the qubits are operated. 
Following the standard technique of adiabatic elimina- 
tion [l4| , we can eliminate the bosonic mode of the large 
JJ and obtain an effective interaction Hamiltonian be- 
tween the two qubits: X\2<Jx^ &x\ with the coupling 
strength X12 = — 2310920 /fko p . Using the expression of 
5io i 520: an d u) p , one can easily confirm that this inter- 
qubit coupling is the same as that in Eq. (J2j. Here the 
large JJ serves as the data bus to virtually mediate the 
interaction between the two qubits. Therefore, the clas- 
sical and quantum treatment to the large JJ are equiva- 
lent to each other. Generalizing the two-qubit system to 
the multi-qubit system, the effective inter-qubit coupling 
term reads J2i> 3 \XijOx a x with Xij = -^9io9jo/^p- 
Because g i0 = - (Eji/2q) sin (7r$ e /$ ), the coupling Xij 
is tunable by changing the static magnetic field <& e ap- 
plied to the loop. We should point out that if the dc 
magnetic flux <E> e is replaced by an ac variable-frequency 
magnetic flux $ e (i), then the qubit can be selectively cou- 
pled to the data-bus by a well-chosen frequency-matching 



condition between the qubit, data bus, and the variable- 
frequency magnetic flux. 



V. CONCLUSIONS 

In summary, we have studied a variable-frequency- 
control approach in charge-qubit circuits: the switchable 
coupling between the two charge qubits can be imple- 
mented by changing the frequency of the externally ap- 
plied magnetic flux. Single-qubit operations can also be 
addressed and operated selectively. The charge qubits are 
chosen to work at their optimal points, so the effect of the 
noise, resulted from uncontrollable charge fluctuations, 
on the charge qubits is much suppressed. Moreover, the 
effects of the microwave field on the supercurrent of the 
two qubits are discussed. It is found that the supercur- 
rent of the qubits significantly depends on the states of 
the microwave field. We also discuss the quantum treat- 
ment of the large JJ and find that both quantum and 
classical treatments are equivalent to each other. If the 
two-qubit circuit is generalized to many qubits, the in- 
teraction Ox&x can also be achieved. 
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